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Quantum systems unfold diversified correlations which have no classical counterparts. These 
quantum correlations have various different facets. Quantum entanglement, as the most well known 
measure of quantum correlations, plays essential roles in quantum information processing. However, 
it has recently been pointed out that quantum entanglement cannot describe all the nonclassicality 
in the correlations. Thus the study of quantum correlations in separable states attracts widely 
attentions. Herein, we experimentally investigate the quantum correlations of separable thermal 
states in terms of quantum discord. The sudden change of quantum discord is observed, which 
captures ambiguously the critical point associated with the behavior of Hamiltonian. Our results 
display the potential applications of quantum correlations in studying the fundamental properties 
of quantum system, such as quantum criticality of non-zero temperature. 
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Nowadays the significance of quantum correlations 
goes well beyond quantum entanglement, which has been 
widely investigated and believed to be the key resource 
of quantum information processing [ ]. However, quan- 
tum correlations possess other facets for which the quan- 
tum entanglement does not provide a complete charac- 
terization: there exist quantum nonlocality without en- 
tanglement [ ]. Moreover, quantum entanglement is not 
a necessary for quantum computation: we presented the 
experiment to perform a quantum algorithm for parity 
problem without using entanglement as early as 2001 [ ], 
and recently a deterministic quantum computation with 
one qubit (DQCl [ ]) is realized experimentally [ ] on 
the mixed separable states. In this context, it is actu- 
ally the type of quantum correlation known as quantum 
discord [ , ], rather than quantum entanglement, that 
provides the enhancement for the computation [ ]. The 
significance of quantum discord in quantum communica- 
tion has been studied in the cases of local broadcasting 
[9] and state merging [10]. 

Quantum correlations also provide a powerful frame- 
work for the understanding of complex quantum systems. 
The first step towards this investigation is due to Os- 
borne and Nielsen [11], who studied the relation between 
the entanglement and the quantum phase transition in 
XY model. Recently, it is shown that quantum discord, 
in contrast to entanglement, captures the critical points 
associated with quantum phase transitions for XX Z and 
XY model, even at finite temperature or in an external 
magnetic field [12-16]. Besides, quantum correlations are 
of great importance in quantum thermodynamics [17, 18], 
and even in photosynthesis [19]. Thus the study of the 
quantum correlations will be of both fundamental and 
practical significance. 

In this Letter, we report an experiment to character- 
ize the quantum correlations in separable thermal states. 



Quantum discord is used as the quantifier of the quan- 
tumness in the correlation in solids. A series of sepa- 
rable thermal states are generated, which can be taken 
as the thermal-equilibrium states of a two-qubit XX Z 
Heisenberg model at a finite temperature. By tuning 
the state parameter that corresponds to the anisotropic 
coupling constant in the XXZ model, we observe the 
sudden changes of quantum discord. The sudden change 
corresponds exactly to the energy-level crossing of the 
XXZ Hamiltonian, and also indicates the critical point 
on which the ground state of the Hamiltonian transforms 
from product state to entangled one, or vice versa. When 
the number of the qubit on the XXZ chain tends to in- 
finity (i.e., thermodynamic limit), the sudden changes of 
discord spotlight the critical points associated with quan- 
tum phase transitions [ ]. Thus our experiment opens 
the possibility of experimentally studying the fundamen- 
tal properties of quantum system from the viewpoint of 
quantum correlation, in particular quantum discord. 

The thermal states generated in our experiment take 
the Bell-diagonal form, that is. 
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where a!^^^^ {i = x, ?/, z) 
(second) qubit. From | 
is given by 



are the Pauli matrices of the first 
] , the quantum discord of pbd 
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where c = max{|c^|, Ic^^l, |c^|}, S{p) = -Tr(plog2p) is 
the von Neumann entropy, and the function h{x) is de- 
fined as h{x) = -^^ log2 ^-^ - ^^ log2 ^-^. 

The Bell-diagonal states of the form (1) can be taken 
as the thermal-equilibrium states of a two-qubit XXZ 



Heisenberg chain with the Hamiltonian given by 
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where J is the couphng constant and A is anisotropy 
parameter. In fact, the density matrix for this model 
at thermal equilibrium at temperature T is given by the 
canonical ensemble p = e~^^/Z, where /3 = l/ksT^ and 
Z = Tr(e"^^) is the partition function. When T = 
0, the state p is the ground state of the Hamiltonian 
H. When the temperature T is appreciably larger than 
the maximum splitting of H in Eq. (3), the thermal- 
equilibrium state can be approximated as 
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The entanglement of formation (EoF) [21] and the quan- 
tum discord of the state p can be easily obtained from 
Ref. [22] and Eq. (2) respectively. We depict the results 
with respect to the anisotropic parameter A in Fig. 1, 
where the units have been chosen such that J = 1 and 
kB is unity [ ]. It is shown that with the temperature 
increasing, the EoF becomes zero, while the quantum dis- 
cord always remains positive [Fig. 1(b)]. Also, at finite 
temperature, the quantum discord changes suddenly at 
A = ±1 which corresponds exactly to the energy- level 
crossing, but it is not the case for EoF: The behavior of 
EoF is illustrated by smooth curves (grey lines in Fig. 

1(b)). 
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FIG. 1. (color online) (a) the energy- level of the Hamiltonian 
given by Eq. (3). (b) the EoF (grey lines) and quantum 
discord (red lines) of the state p given by Eq. (4). When the 
temperature increase from OK to 2K, the EoF vanishes and 
the discord remains positive. The sudden changes of discord 
correspond to the level crossing shown in (a). 

To explore the quantum correlations in the separable 
states, we choose phosphorous donors in silicon (P:Si) 



material [24, 25] with P concentration about 1 x lO^cm"^ 
as a benchmark system. The system consists of an elec- 
tron spin /S = 1/2 and a nuclear spin 7 = 1/2. It can be 
described by an isotropic spin Hamiltonian: 
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where uOe = gPeBo/h and ujj = giPjBo/h character- 
ize the Zeeman interaction for the electron and nuclear 
spins, a = 117 MHz is the isotropic hyperfine interac- 
tion strength and S ( I ) is the electron (nuclear) spin 
operator. The energy diagram of this system is plot- 
ted in Fig. 2(a), where the two electron-paramagnetic- 
resonance (EPR) transitions and two NMR transitions 
have been labeled by MWl, MW2 and RFl and RF2 
respectively. The frequencies of these transitions, which 
are determined by field-swept echo detection (FSED) and 
standard Davies electron-nuclear double resonance (EN- 
DOR) experiments, are MWl = 9.701 GHz, MW2 = 
9.818 GHz, RFl = 52.383 MHz and RF2 = 65.181 MHz. 
The pulse length of the EPR (NMR) tt pulse is 60 ns 
(10 ps) which has been derived from the EPR (NMR) 
Rabi nutation experiment. 



The experiment has been carried out at Te: 



K on 



a homebuilt pulsed electron spin resonance (ESR) spec- 
trometer, which provides the access to controlling both 
the electron and nuclear spins with fiexible microwave 
(MW) and radio frequency (RF) pulses [26]. We first 
measure the relaxation times of both electron and nu- 
clear spins (see Fig. 2(b)). For the electron spin, the 
transverse relaxation time is T2e = 120 ps (green line in 
Fig. 2(b)) and the electron population relaxation time 
is Tie = 5.6 ms (red line in Fig. 2(b)). For the nuclear 
spin, the dephasing time is determined as T2*^ = 24.3 ps 
by a nuclear spin free induction decay (FID) experiment 
(blue line in Fig. 2(b)). We also measure the decay of 
the p23 under the environment with a decay constant 
tc = 200 ns (plotted in the inset of Fig. 2(b)). The re- 
laxation times are similar to those reported in Ref. [24]. 
The nuclear population relaxation time Tin is estimated 
to be 250 times of Tie according to the literatures [27]. 
So the waiting time between each experiment has been 
set to 10 s. 

Fig. 2(c) shows the pulse sequence applied in our ex- 
periment. The first MW pulse is to fiip the electron spin 
with an angle while Iz = —1/2. There is a waiting 
time (ti = 1 ps ^ T|g) following so that the off-diagonal 
elements of the density matrix decay off. The first 90° 
RF1(RF2) pulse is to equalize the diagonal elements pu 
and P22 (P33 and ^44). Then we let the off-diagonal ele- 
ments, which are generated by the first RF pulse, decay 
off by another waiting time r2 = 200 ps which is long 
enough compared with T^^. The following 90° RF2(RF1) 
pulse (90°) and the 180° MW1(MW2) pulses are to gen- 
erate the nonzero off-diagonal element, p23- After waiting 
time Ts which is to introduce a decay factor A(r3), the fi- 
nal state has been prepared to Bell diagonal states given 
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FIG. 2. (color online) (a) Energy level diagram of the P:Si sys- 
tem. There are four Zeeman product states which are labeled 
states 1—4, where the left and the right t (i) stands for the 
zbl/2 states of electron and nuclear spins, respectively. EPR 
and NMR transitions are indicated by two-way arrows, (b) 
Experimental results of the relaxation time measurements. 
The red line are collected from the echo-detected inversion 
recovery experiment (the related time axis is the red bottom 
one), which reveals Tie. The green line is the electron spin 
echo decay through the Hahn echo sequence, which gives T^e- 
The blue curve is the result of the nuclear spin FID with pulse 
sequence (7r2,4 — 7r/2i,2 — t — 7r/2i,2— echo, r a variable), which 
is used to calculate T2n • Note that the time axis for the green 
and blue curves is the black top one. The inset shows the de- 
cay of P23 • The experimental data (black circles) were fitted 
with a damped exponential function (grey line) of the form 

y = yo -\- ^e *c ^ . (c) Diagram of the experimental pulse 
sequence. We use upper (lower) pulse sequence to generate 
the state p, in the case of c^ > {cz < 0). 62,4 (180i,3) 
is a MW1(MW2) pulse to rotate the electron spin by (tt). 
90i,2(3,4) stands for a 7r/2 RF1(RF2) pulse. The tomography 
procession is composed of readout pulse sequence and electron 
spin echo detection. 



by Eq. (1), where Cx = Cy = — e(l — cosO)X{rs)^ Cz = 
±2e(l + cos6>) and e = g(3eBo/8kBTe^p = 7.35 x IQ-^. 
Here X{rs) = exp(— rs/tc) which has been determined in 
Fig. 2(b). 

Tomography technique [28] is then used to reconstruct 
the density matrix p. By observing both the in-phase 
and quadrature components of the electron spin echo 
(the block labeled by electron spin echo detection in Fig. 
2(c)), measurement in Sx and Sy bases is achieved. For 
the diagonal elements pu^ i=i-4, the readout pulse se- 
quence is a MW1(2) (RFl (2)) pulse (green block in Fig. 
2(c)) with variable duration. Note that an additional 
waiting time r4 = 200 ps (not shown in Fig. 2(c)), which 
is long enough to let dephasing effect eliminate the pos- 
sible off-diagonal elements, is inserted before the readout 
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FIG. 3. (color online) (a) The readout of the p23- l-\-x—l-x 
(red circles) and I+y—I-^ (green triangles) are plotted. Sine 
functions (dashed lines) are fitted to the experimental data. 
The X axis is the rotation angle (j) and the y axis is the ampli- 
tude of the echo in unit of e. The results give both real and 
imaginary parts of p23. (b)The reconstructed deviation den- 



sity matrix dp = p- 
and Cz = 0.0008. 



^1 in unit of e, when Cx 



-0.0044 



pulse. Then this selective electron (nuclear) spin Rabi nu- 
tation corresponds to the measurement of SzI'^'^{S^'^Iz) 
where S'(/)"'^ = (1 ± cfz)/2. To readout the off-diagonal 
elements /)i2, P34, P13 and p24, T4 is set to zero. Compared 
with the previous nutation experimental data, these off- 
diagonal elements were calculated to be almost zero in 
our experiment. 

To measure p23, the readout pulse sequence is com- 
posed of a microwave pulse, 1802,4, and a radio fre- 
quency pulse, 03,4, with variable rotation angle (j). We 
denote the spectrum of Rabi nutation, when the phase 
of 1802,4 is ±x, ±?/, as I±x,±y{(l))' After some calcula- 
tions, we find that /+a^((/)) — I-xW = — ^e(p23) sin(0) 
and /+2y((/)) - I-y{(l)) = /m(p23) sin(0) [Fig. 3(a)]. For 
the measurement of pi4, we use ^1^2 instead of (j)s,4: in the 
readout pulse sequence of P23. To normalize the results, 
we compare all the the amplitudes of the Rabi nutations 
in our tomography procession to the amplitude of elec- 
tron spin Rabi nutation which can be taken as 2e. Once 
the above steps are completed, the experimental density 
matrix can be fully reconstructed. Fig. 3(b) shows the 
result of the tomography of one Bell diagonal state whose 
Cx = Cy = -0.0044 and c^ = 0.0008 (Note that 6p is 
the result the reconstructed density matrix subtracted 
by 1/4). The entanglement of this quantum state is zero, 
thus the state is separable. However, the quantum dis- 



cord is 1.45 x 10~^. This shows that nonzero quantum 
correlations have been generated in a separable quantum 
state. 
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meaningful step to observe quantum criticality at finite 
temperatures via quantum correlations. 
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FIGl 4. (color online) Ihe behavior of quantum correlations 
when A varies. The x axis is the tuning parameter A and the 
y axis is quantum discord. Experimental data of quantum 
discord (red circles) , which are derived from the reconstructed 
densfty matrices, agree with the theoretical prediction usmg 



Eq. (2) (solid line). The deviation between the experimental 



data! and the theoretical prediction is due to the impertection 



of the pulses. 



T he values of quantum diacord (red circlca in Fig. 4(b)) 
are numerically calculated uaing Eq. — {^. — The cxpcri 
mental data agree with the theoretical prediction (line in 
the Fig. 4(b)). It is clear that quantum discord experi- 
enc es a sudden change at A = ±1 where the energv level 
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cros sing occurs. It is worthv of pointing out that the en- 



tan glement, which is characterized b}^ EoF, is measured 



to bb zero in our experiments — Thus, we have presented 



a case where entanglement fails to capture energy level 
crossing while temperature is of a finite value. Note that 
there are theoretical works which conclude that quantum 
discord can be utilized to capture the quantum phase 
transition even at finite temperatiires[14, 16, ?,9i\, The 



succless in highlighting the sudden change of the ground 



state (A = — 1) in the two-qubit system employed in our 
expfi'lmenlri IndicaLori LhaL quanLum coiTelaLion could bo 



used Lo obrit^i'Vt^ Llit^ quaiiLum ci'lLlcaULy aL flnlLe Lemper- 
atures. 

In summary, we have experimentally observed quan- 
tum correlations in a series of separable states in solids. 
Furthermore, the capability of using quantum correla- 
tions to reveal the intrinsic change of the physical system 
is rqvealed. Two-qubiL XXZ Heisenberg model has been 



12 



shmj@ustc.edu.cn 
djf@ustc.edu.cn 
R. Horodecki, P. 



Horodecki, M. Horodecki, 



and 



K. Horodecki, Rev. Mod. Phys., 81, 865 (2009) 



C. H. Bennett, D. P. DiVincenzo, C. A. Fuchs, T. Mor, 
E. Rains. P. W. Shor. J. A. Smo lin. and W. K. Wootters, 



Phys. Rev. A, 59, 1070 (1999) 



J. Du, M. Shi, X. Zhou, Y. Fan. B. Ye . R. Han, and 



J. Wu, [Phys. Rev. A, 64, 042306 (2001). 

E. Kni ll and R. Laflamme, [Phys. Rev. Lett., 81, 5672| 



(1998) 



B. P. L anyon, M. Barbieri, M. P. Almeida, and A. G. 

White, |Phys. Rev. Lett., 101. 200501 r2008l|. 

H. Ollivier and W. Zurek, [Physical Review Letters, 88, 
017901 (2002)| . 



L. Henderson and V. Vedral, [J. Phys. A: Math. Gen., 34, 
6899 (2001)| . 



A. Datta, A. Shaji, and C. M. Caves, Phys. Rev. Lett.,| 



100, 050502 (2008)| . 



M. Piani, P. Horodecki, and R. Horodecki, Phys. Rev 



Lett., 100, 090502 (2008) 



D. Cavalcanti, L. A olita, S. Boixo, K. Modi, M. Pi- 
ani, and A. Winter, [Phys. Rev. A, 83, 032324 f2011)| ; 



V. Madhok and A. Datta, N6zc?., 83, 032323 (201ll 



T. J. Osborne and M. A. Nielsen, phys. Rev. A, 66, 
032110 (2002] . 



R. Dil lenschneider, [Physical Review B, 78, 224413| 



(2008) 



M. S. Sarandy [Phys. Rev. A, 80, 022108 (2009) 



T. Werlang, C. Trippe, G. A. P. Ribeir o, and G. Rigolin, 
Phys. Rev. Lett., 105, 095702 (2010) . 



J. Maziero, H. C. Guzman, L. C. Celeri, M. S. Sar andy. 



and R. M. Serra, |Phys. Rev. A, 82, 012106 (2010) 

T. Werlang, G. A. P. R ibeiro, and G. Rigolin, |Phys. Rev 



A, 83, 062334 (2011 



W. H. Zurek, [Phys. Rev. A, 67, 012320 (2003)[ 

K. Maruvani a, F. Nori, and V. Vedral, [Rev. Mod. Phys.,| 



81, 1 (2009) 



M. Sarovar, A. Ishizaki, G. R. Flem ing, and K. B. Wha- 
ley, [Nature Physics, 6, 462 (2010) . 



S. Luo, [Phys. Rev. A, 77, 042303 (2008) 



C. H. Be nnett, H. J. Bernstein, S. Popes cu, and B. Schu- 



macher, Phys. Rev. A, 53, 2046 (1996) 



W. Wootters, Physical Review Letters, 80, 2245 (1998). 
X. Pen g, J. Du, and D. Suter, [Phys. Rev. A, 71, 012307 



takeii as aii example, and the abrupt change of ILs ground 
state has been unambiguously spotlighted by quantum 
correlations. Our experiment may serve as a preliminary 



(2005)[; X. PenR, S. Wu, J. Li, P. S uter, and J. Du 



Phys. Rev. Lett., 105, 240405 (2010) 



[24] Y. Wang, X. Rong, P. F eng, W. Xu, B. Chong, J.-H 



Su, J. Gong, and J. Du, Phys. Rev. Lett., 106, 040501 



[25] 



(2011). 



G. Briggs, J. Ager, and S. Lyon, Journal of Physics 



A. M. [Tyryshkin, S. A. Lyon, A. V. Astashkin and A. M. 



Raitsimring, Phys. Rey. B, 68, 193207 (2003) 



[26] J. Du, M. Shi, J. Wu X. Zhou, and R. Han, 
|A, 63, 042302 f200ll|. 

[27]'- - --■ - -- 



Condensed Matter, 18, S783 (2006) 



Phys. Rey. 



[28] L. Vandersypen and I. Chuang, Reyiews of Modern 
Physics, 76, 1037 (2004) 



[29] T. We rlang and G. Rigohn, Phys. Rey. A, 81, 044101 



A. Tyryshkin, J. Morton, S. Benjamin, A. Ardayan, 



(2010) , 



